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Abstract 

In this paper we give fully analytical description of the dynamics of a collection of 
7V-Frenkel excitons in high density regime dispersively coupled to a single mode cavity 
field, in the presence of both exciton and cavity-field dissipations. By using excitonic 
operators as g-deformed bosonic operators for the system, we solve analytically the 
Liouville equation for the density operator at zero temperature and investigate the 
influence of the number of excitons and the effect of both dissipations on dynami- 
cal behavior of the system. We use the solution of master equation to explore the 
dissipative dynamics of non-classical properties such as, molecule-field entanglement, 
quadrature squeezing of the field, and molecular dipole squeezing. We find that the 
non-classical properties are strongly affected by the number of excitons and also by 
the existence of both dissipations. 

1 Introduction 

In quantum theory and cavity-quantum electrodynamics (Cavity-QED) of semiconductors 
where cavity is chosen as a bulk or confined systems such as quantum wells and quantum 
dots, we work in the low density regime of excitons, and treat these quasiparticles as ideal 
bosons to solve wide number of problems in quantum optics, cavity-QED and quantum 
information [H |2]. Typically, these studies have been done under the condition that 
the influence of the environment is not taken into account. The environment which is 
represented by a thermal reservoir always exist, and affects the system under consideration. 
No matter how weak the coupling to such an environment, the evolution of quantum 
subsystems is eventually affected by non-unitary features such as decoherence, dissipation 
and heating. The initial information irreversibly leaks out the system into the very large 
number of uncontrollable degrees of freedom of the environment. Initially prepared pure 
states are typically corrupted on extremely short time scales due to quantum coherence 
loss that turns them into mixed states [3]. The disspersive effects caused by the energy 
exchange between the system and environment have been studied in the Jaynes-Cummings 
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model (JCM) [4j, both analytically [5j and numerically "6]. In the last few years the JCM 
with phase damping, as applied to decoherence and entanglement [7] , dissipative dynamics 
of JCM model [8], and nonlinear quantum dissipation effects on dynamical properties of 
the f-deformed JCM model [9] have been studied. In addition to the exact solvability of 
JCM within the rotating wave approximation, one of the most interesting aspects of its 
dynamics is the entanglement between atom and the field. Entanglement as a physical 
resource has been used in quantum information science such as quantum teleportation 
[To] , super dense coding \TT\ and quantum cryptography [T2] . It has been shown that due 
to the influence of the field dissipation in JCM under dispersive approximation that the 
amplitude of the entanglement between the field and the cavity-field decreases with the 
time and at least completely suppresses [8l[9l[13]. Sub-Poissonian statistics and quadrature 
squeezing of the cavity-field are two remarkable non-classical effects have also been studied 
in different versions of generalized JCMs [9, 14J. 

Fascinating electronic and optical properties of spatially confined nanostructures like 
quantum wells and great potentially of such structures in semiconductor cavity-QED ap- 
plication has motivated permanent extension of their study. Among a variety of new 
results in this field, it is important to find out the influence of dissipation on non-classical 
properties, so in our model, we replace one atom in JCM with A^-Frenkel excitons. How- 
ever, if we want to replace two-level atoms in the cavity with high density approximation 
of excitons in semiconductor nanostructures, we can't use the standard Dicke model be- 
cause the excitons are not more ideal bosons. One way to deal with the problem is to 
replace ideal bosons with g-deformed ones in the Dicke model [15] . We also consider the 
influence of dissipation between excitons and their reservoir, visualized as a large number 
of harmonic oscillators, and also the influence of cavity-field damping due to its coupling 
to the environment, on the dynamics of the system. For this purpose, we use the Liouville 
equation for the density operator to study the dynamical behavior of a dissipative system 
composed of iV-Frenkel excitons (iV-two level molecules) interacting with a single mode 
cavity-field and we show the influence of dissipations and the number of molecules on 
quantum optical properties of system such as, entanglement between the cavity-field and 
molecules, quadrature field squeezing and also molecular dipole squeezing. 

This paper is organized as follows. In Sec. 2 we first introduce the Hamiltonian of 
the system without considering any dissipation in the presence of excitonic operators as q- 
deformed bosons. In Sec. 3, we give an analytical solution for the master equation to derive 
the total density operator in the presence of both dissipations. In Sec. 4 we employ the 
analytic results obtained in section. 3 to investigate the influence of number of molecules 
and the effect of both reservoirs on the dynamical properties of molecules and cavity-field. 
Finally, we summarize our conclusion in Sec. 5. 

2 The Hamiltonian of the system 

We consider our system as N identical two-level molecules of splitting uJeg interacting 
with a single-mode quantized cavity- field of frequency cvq. We assume that all molecules 
have equivalent mode position, so they interact with the cavity-field by the same coupling 
constant g. In the case of relatively high density of molecules in the excited state the 
Hamiltonian of the system can be written in terms of the exciton operators bq,bq (as 
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g-deformed bosonic operators) [15] 

H = blbqcveg + atJo + VNg (^d^q + bla^ , h=l, (2.1) 

where a and are the annihilation and creation operators of the quantized cavity-field 
respectively, and the g-deformed bosonic operators bq and bl satisfying the g-deformed 
commutation relation 

KblU = bqi\ - Qblbq = 1, 9 = 1-^, (2.2) 

are defined as bq = b^b^j \fN , b\ = bgb^ j V^, in which be, bg {bl, bg) are the annihilation 

and creation (creation and annihilation) operators of one molecule in the excited state 
and the ground state, respectively pS]- Following Peixoto et all, [8j, in the large detuning 
approximation (dispersive limit), that is 



— > Jriph + 1, 6 = u;eg- Wo, (2.3) 
9 

for any relevant photon number Uph , we obtain the following effective interaction Hamil- 
tonian 

H^jj- = ujAm,n (n,m| (d^d -|- 1) -|- d^d |n — 1, m -|- 1) (n — 1, m -|- 1| | , (2.4) 

where w6 have defined. -A^yj — 

n{m + 1), and u> = g'^/5. The state \n,m) indicates that 
n molecules (m = N — n molecules) are in the excited state (ground state). Furthermore, 
we have used 



bq \ n, m) 



b^ \n, m) 

The effective Hamiltonian ()2.4p which does not cause any transition in the system creates 
an entanglement between the molecular and the field states. 



1 t' ft I Jnim -|- 1) 

— 6e&J |n, m) = — |n - 1, m + 1) , 



b\b. 



N 



>g \n,m) 



m 



\n + 1,171 — 1) 



(2.5) 
(2.6) 



3 The master equation and its analytical solution 

We assume that our system is coupled with two types of reservoirs, one of them is coupled 
to the cavity-field and the other is coupled to molecules. The time evolution of the 
collection of N two-level molecules interacting dispersively with a single mode cavity-field 
in the interaction picture can be described by the following master equation 

^P{t) = -i [^iff-iPi^)] + D molecules P{i) + bfieldP{t), h = 1, (3.7) 
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where p{t) is the density operator of the system, and belongs to the set D{'Hm 'Hp) of 
the trace class operators acting in the space corresponding to the direct product of the 
two Hilbert spaces Tiu and Hp of the molecules and the field, -fffjj is given by (j2.4p . 
and the damping of the cavity-field and molecular subsystems are phenomenologically 
represented by the superoperators Dji^id and Dmoiecuies- By using thermal reservoir at 
zero temperature for the field and one atom |16p i7| and the same one for molecules, using q- 
deformed bosonic operators instead of atomic operators, we can write these superoperators 
as 

D field' = k{2a ■ — a) a ■ — ■ a^a), (3-8) 

Dmoiecuies- = k' {2\ ■ b\ - b\hq h\bq), (3.9) 

where k and k' are the corresponding dissipation constants. These superoperators are 
linear combinations of bosonic superoperators, and form a finite Lie algebra under com- 
mutation. The bosonic superoperators represent the action of creation and annihilation 
operators of the harmonic oscillator on an operator O : 

(a.)0 = a^O = ad, (at.)O = (a^)tO = a^O, , 
(.a)0 = a'O = da, (.at)0 = {a'^O = Oat. ^> 

Similarly, the action of g-deformed creation and annihilation operators bq and bq on an 
operator O can be defined as 



{bq.)d = bid = bqd, (11)6 = {bD^d = Pqd, 
{.hq)d = bp = dbq, {:bl)d = (6pto = dbl 



(3.11) 



To solve the master equation (|3.7p . we rewrite it in the basis of molecular states \n,m) 
(note that there is just one mode of the cavity-field in the cavity which can interact with 
one molecule) 

Pn,n Pn,n— 1 \ / -^n,n -^n,n~l 



Pn~l,n Pn—l,n—l J \ -^n—l,n 



'-n— l,n— 1 



(3.12) 



where p{t) = 4zp{t), and we have defined 



dt' 

Pn,n{t) = {n,m\ p{t)\n,m) , (3.13a) 

Pn,n^iit) = {n, m\ p{t) \n-l,m + l) , (3.13b) 

Pn-i,nit) = {n-l,m+l\ p{t) \n, m) , (3.13c) 

Pn-i,n-i{t) = {n-l,m + l\ p{t) \n-l,m + l) , (3.13d) 



Furthermore 

Xn,n = [iAmMP " ^) " 25n,m^' + K'^J 'M- P)} /5„,„(t), (3.14a) 
Xn,n-1 = {-iAm,nOj{M + P + 1) - Bl^^k' + k{2J - M - P) } (t) , (3.14b) 

= [iAmM^ + P + 1) - Bl^k' + k{2J -M-P)} p„_i,„(t), (3.14c) 

Xn-l,n-l = {iAm,nUj{^ " -P) + H^J -M-P)} Pn-l,n-l{t) + 2Bl^k' Pn,n{t) , 

(3.14d) 
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are operators that act in TCp (remember that n shows the number of excited molecules). In 
Eqs. (|3lil) . we have defined Bn,m = \/n{m + 1) j \/]V = ^ 

m,n /\/^) and also M = a}a- 

, P = -a^a , J = a ■ a) , which satisfy the following commutation relations 

[J, M] = J, [J, P] = J, [M, P] = 0. (3.15) 

In order to find the total density operator we need to solve the four Liouvillian equations 
associated with each matrix element. For any initial state, the total density operator 
can be evaluated from the solution of Eq (j3.12p . We assume that the initial state of the 
interacting system is given by 



= —i= m) + |n - 1, m + 1)} (g) \a) , 
v2 



(3.16) 



which means that at t = 0, molecules are in the superposition of two states |n, m) and 
|n — l,m + 1), and the cavity-field is prepared in the Glauber coherent state. After some 
calculation we evaluate the matrix elements of the field density operator, using initial 
condition (|3.16|) . as 



Pn,n(i) = 
Pn,n-l(i) — t; 



(3.17a) 
(3.17b) 
(3.17c) 



Pn-l,n-l(i) - 2 



1 



(ae~'^*e*^'"'"'^*)" (Q,*g-'=*g-*^m,n^*)"' 



n,n' 

X |^g-2Bn,mfc'tg2i^m,„^t(n'-n) _ ^ 

where the definitions 

T{t) = -|a|2(l-e-2'=*) 



n 



){n'\ 



2R2 U 



-2Bl k' + 2iAm,nOj{n' - n) 



(3.17d) 



I |2 7 



e ^''*(A;cos(2Am,„wt) - A^.nW sin(2Am,nWt)) - fc) , (3.18) 



(3.19) 



have been used. As a result, the total density operator of the interacting system can be 
written as 

Ptotai{t) = Pn,n{t) ® |n, m) (n, m\ + Pn,n-i{t) ® |n, m) (n - 1, m + 1| 

-|-pn_i,n(t) (g) |n - 1, m + 1) (n, m| + /5n_i,„_i(t) (g |n - 1, m + 1) (n - 1, m + 1| . 

(3.20) 
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By taking the trace of the total density operator (j3.20p on the molecular (field) vari- 
ables, we get the reduced field (molecular) density operator: Pf(t) = Tr molecule {ptotai{i)) 
( Pm{t) = Trfi^id{ptotai{t)))- In the next two sections we are going to use the solution 
given by ()3.20p to investigate the dynamical properties of the system. 



4 Dynamical properties of the model 

4.1 Linear entropies and molecules-field entanglement 

It is well-known that if the cavity-field and molecules are initially prepared in a pure state, 
then at i > the molecules-field system evolves into an entangled state. In this entangled 
state the field and molecules are separately in mixed states. The stability of a pure state 
may be understood as the process where quantum coherence of the state is preserved along 
its time evolution. In this sense we say that an initial pure quantum state, described by 
the density operator p{t) is stable if Trfp'it) = 1, for all times. One way to measure the 
stability of an initial pure state is to use the linear entropy [18j 

s = l- Trp^{t). (4.21) 

The time evolution of the molecule(field) entropy reflects the time evolution of the degree 
of entanglement between the molecules and the field. The higher the entropy is, the greater 
the entanglement between the molecule and the field becomes. By using Eqs. (j3.20p and 
(j4.2ip the linear entropy of the total system under consideration is obtained as follows 

sm-}{t) = l-Tr {pl,Jt)) = 1(1 + e-^^"-^-'* + 2e2re-2^i™'='* 

n,n' 

+4Am,nUjBlmk'{n - n')e~2^'-'='* sin(2^„,„u;(n - n')t))) ^^'^^^ 

n,n' 

X ^e-^Bl^k't _ 2^-2Bl^k't cos{2Am,nio{n - n')t) + l}). 
Furthermore, the linear entropy of the cavity- field is given by 

Sf{t) = 1 - Trf (p}{t)^ = 1 + e"^^".™'''* + 2e"^'^".'"'''*e2|"l''^"''"('=°"(2A™."'^*)-i) 

X je-^^".™'''* - 26"^-^".'"'''* cos(2ylm,„a;(n' - n)t) + l| 

n,n' 

x{{-{2Blmk'? (e-2^--^'*cos(2^^,„u;(n' - n)t) - l) 
+4^^ ^uj'^Bl mk'{n' - n)e-2-S'.™*^'* sin(2A„,„w(n' - n)t))). 



n,n' 
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and the molecular coherence loss will be measured by 

Sm{t) = l-Trm {pl{t)) 

(4.24) 

According to Eqs. (|4.22[ )-( l4.24| l. we plot the time evolution of the linear entropy of the 
total system, the cavity-field and the molecular subsystem as functions of scaled time wt, 
in Figs. 1. 

First, we discuss the coherence loss of the field. For all diagrams in Figs. 1, the 
field is initially prepared in the coherence state and Sf{t) = 0, but in a short stage 
Sf{t) increases and the coherence of the field loses (of course the equilibrium state of 
the filed correspond to the vacuum, for which the linear entropy of the field is zero). 
The linear entropy of the field Sf shows local maxima and minima, corresponding to the 
entanglement and disentanglement between the field and excitons. In the presence of 
molecular dissipation (Figs, la, lb), the field linear entropy still exhibits local maxima 
and minima, but minimum of disentanglement is not taking place at = and the 
field is not in a pure state. In the absence of molecular dissipation, in figure. Ic, when 
disentanglement happens periodically, the field is in a pure state {sp = 0). 

To verify the role of two types of dissipations on Sm and Sm-f, we discuss the coherence 
loss of the molecules and of the system. We find that the linear entropy of the system 
and of molecules exhibit periodic behavior, and these periods coincide with the related 
periods of sj. Note that the field and molecules are not in pure state. In the course of 
time evolution, time, both Smit) and Sm-fit), in the presence of molecular dissipation 
(Figs, la, lb), find the maximum values (corresponding to maximum coherence loss ), 
and then these linear entropies decrease smoothly to find their asymptotic values at zero 
(the equilibrium state of molecules corresponds to the state that all molecules are in the 
ground state and the linear enrtropies are zero). Furthermore, in the absence of molecular 
dissipation, these two entropies tend to an asymptotic non zero value (Figs. Ic, Id). We 
also show the linear entropies for one molecule in the absence of molecular dissipation 
(Fig. Id), and find the same results as the usual JCM given in reference [8j. 

4.2 Quadrature squeezing of the cavity-field 

In the past two decades, there has been major studies focused on the fluctuations in the 
quadrature amplitude of the electromagnetic field to produce squeezed light. This light is 
indicated by having less noise in one field quadrature than vacuum state with an excess of 
noise in the conjugate quadrature such that the product of canonically conjugate variances 
must satisfy the uncertainty relation. In the studies of quantum optics theory, this light 
occupies a wide area because of the various applications, e.g., in optical communication 
networks [IQj, in interferometric techniques [20j, and in optical waveguide tap [21j. Fur- 
thermore, investigation of the squeezing properties of the radiation field is a central topics 
in quantum optics and noise squeezing can be measured by means of homodyne detection 

In order to investigate the quadrature squeezing of the model under consideration, 
we introduce the two slowly varying Hermition operators Xia{t) and X2a{t) defined by. 
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respectively, 

^ia{t) = \ (ae^'^°* + ate-^'^°*) , X2a{t) = ^ (^ae'^°' - a^e-"^°') . (4.25) 

A state of the field said to be squeezed when one of the quadrature components Xia{t) 
and X2a{t) satisfies the relation 

AXUt)y^<l, (i = lor2). (4.26) 

The degree of squeezing can be measured by the squeezing parameter Si{i = 1,2) defined 

by 

Si{t) = A(^(^AXia(t)y^-l. (4.27) 

Therefore the condition for squeezing in the quadrature component can be simply written 
as Si{t) < 0. In Figs. 2, we have plotted the squeezing parameter si{t) given by 

s,{t) = 2e-2|«l'^-''=* E i^^^^n + 2{e-2^'.-'='* |a|=^ e'^'^* cos(2(a;o - ^m,ncu)i) 

n 

+ |a|2 6-2^=* cos(2(a;o + Am,nUj)t) ^^^,^^^,y,J^2^^^^.,^_^y, 

X {-{2Bl^kr (e-2^"-'='*cos(4A^,„c.i) - ij + S^^^^c^fi^^^fe'e-^^-'-'^'* sin(4A^,„a;i)} 
+ sin(2(u;o + ^m,nC^)t) (^b^ ,^fc0^+4<„c.^(-2)^ 

X [-SAl^^uj'^Bl^k' (e-2«^.'"^''*cos(4A„,,a;0 - l) - {2Bl^kTe-^''lrn'^'' sm{4Am,nCot)}} 

+4{e~2^".™'''*ae"''* cos((tc;o - Am,n^)t) 

+ae-'=* COS((WO + Arn,n(^)t)j^^r^jy)^^U^ 

X {-(2i?2 „^fe')^ (e-2^"-^''*cos(2A„,,,c.t) - l) +4A^,X<mfe'e-2<'"'='*sin(2A^,^^^^ 

+ae-^* sin((^o + ^r»,n^)0 (2i?;i,^fc')^V4A-^,.^^ 

X {-AAl^^u'Bl^k' (e-'^l-''^ cos{2Am,n00t) - l) - (2i?2_^fc')'e-^^"-''* sin(2A^,„a;i)}}' 

(4.28) 

versus the scaled time ut for the corresponding data used in Figs. 1. As it is seen, the 
squeezing parameter si{t) shows fast oscillations and the quadrature component Xia does 
not exhibit squeezing. In the presence of both dissipations (Fig. 2a), squeezing parameter 
si{t) present local maxima and minima due to the field interaction with molecules. In the 
presence of molecular dissipation (Figs. 2a, 2b), the field is no longer in a pure state, except 
initial time and equilibrium state (corresponding to the vacuum), and in the absence of 
molecular dissipation (Fig. 2c), periodically, the field can be find in pure state. 

4.3 MoleculEir dipole squeezing 

In this section we are going to discuss about the dynamical behavior of excitonic ensemble 
property that is molecular dipole squeezing. 
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To analyze the quantum fluctuation dipole variables, we consider dx and dy corre- 
sponding to the dispersive and absorptive components of the amplitude of molecular 
polarization |22j 

<^x(i) = ^(^Je^'"^^* + 6ge^"-*), Oy{t) = ^(ftje"*"-* - 6^6*^-*). (4.29) 

The fluctuations in the component Uiii = x or y) are said to be squeezed if the variance 
in di satisfies the condition (^{Aax{t))'^'^ < | \ {o'z{t))\ {i = x or y). Since (o'iit)) = 1/4 
the condition of dipole squeezing may be written as 

Fiit) = l-i{aiit)f -\{a,{t))\<0 {i = xory) (4.30) 

We find Fy{t), corresponding to the squeezing of dy{t) for the system under consideration, 
as follows 

Fy{t) = 1-4 I lB„,^ere-^'.'"'='*e-l"l''^~'"(^°^(2A„--*)-i) sin(a;e3t + - \a\^ e'^^'* sin(2A^,„u;t)) 

(4.31) 

To investigate the influence of both kinds of dissipations, we plot the time evolution of 
Fy{t) in Figs. 3-6. As it is seen, in the course of time evolution, the function Fy{t) 
shows rapid oscillations which are damped in the presence of dissipations, and under the 
condition that i?n,m 7^ 1, dipole squeezing can occur. In Figs. 3, we show the influence 
of dissipations on Fy{t), in the presence of both dissipations (flgure 3a), in the absence 
of fleld dissipation (flgure 3b), and in the absence of molecular dissipation (flgure 3c). 
We can see that with suppression of dissipations in Figs. 3b, 3c, the rate of damping in 
oscillations decreases, and in the absence of molecular dissipation the oscillations continue 
to narrow band between two positive numbers. In Figs. 4a-4c, we show the influence of 
total number of excitons N on Fy{t). By increasing the number of excitons the squeezing 
can be much stronger and the oscillations damped faster. 

Note that in all diagrams, under the condition Bn,m 7^ !> dipole squeezing appears and 
for any total number of excitons A^, as it is clear in relation ()4.3ip . for larger numbers of 
Bn,m: the effect of squeezing will be stronger. In Figs. 5a-5c, we show the influence of the 
parameter -Bn,m on Fy{t). As we can see, with -Bn,m = 1 there is no squeezing for Fy{t). 
There is a an interesting instance for when molecular dissipation can be suppressed that 
by increasing the number of molecules, Fy{t) shows squeezing all the time, periodically 
(figure. 6a,6b), and by increasing the number of molecules it will be much stronger. 

5 Summary and Conclusions 

In this paper, we have studied theoretically the influence of the dissipations on non-classical 
properties of the system as A^-Frenkel excitons interacting with a single mode cavity-fleld 
in the dispersive approximation. Using g-deformed bosonic operators for excitons, we ob- 
tained an effective interaction Hamiltonian and we used it to flnd the master equation 
containing molecular dissipation and cavity-fleld dissipation. By solving the master equa- 
tion, the total density operator was obtained and some of non-classical properties were 
studied. In continue we summarize our conclusions: 
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First, we discussed the linear entropy to find the coherence loss of the cavity- field Sf(t), 
molecules Sm{t) and of the system Sm-f{t). We found that the linear entropies present 
local maxima and minima corresponding to entanglement and disentanglement due to the 
field interaction with excitons. As the time goes on, in the presence of both dissipations, 
all linear entropies decrease smoothly to find their asymptotic zero values. 

second, we discussed about quadrature squeezing of the field. We found that in the 
presence of dissipations, squeezing parameter si{t) shows fast oscillations with local max- 
ima and minima due to the field interaction with molecules and the quadrature component 
Xia does not exhibit squeezing. 

In the last section we discussed the influences of number of molecules and both dis- 
sipations on the temporal evolution of molecule dipole squeezing Fy{t). We found that 
Fy{t) shows rapid oscillations between positive and negative values that the minus values 
corresponding to the squeezing of ay{t) and Squeezing can be occur when Bn^m 1 (the 
parameter i3n,m has introduced in section 3, is related to the number of molecules which 
are in the excited and ground state), and so there is no squeezing when we have just one 
molecule in the cavity. We also found that by increasing the number of molecules, the 
squeezing can be much stronger and oscillations of Fy{t) damped faster. 
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Figure 1. Time evolution of the linear entropy Sa-f{t) (black solid curve), Sf{t) (gray 
solid curve) and Sa{t) (dashed curve) as functions of scaled time ut, for a = 1, A/" = 10 
and (Fig.la) k = k' = O.OSu;, n = m = 5 (Fig.lb); k = 0, k' = 0.05a; (Fig.lc); k = 0.05a; 

k' = 0, An,m = Bn,m = 1 (Fig. Id). 



Figure 2. The time evolution of si as a function of scaled time ut, for N = 10, 
n = m = 5 a = 1, and for A; = fc' = 0.05a; (Fig.2a); k = 0, k' = 0.05a; (Fig.2b); k = 0.05a;, 
k' = (Fig.2c). 



Figure 3. The time evolution of Fy{t) as a function of scaled time tot, for N = 10, 
m,n = 5, a = 1, and when Bn,m 1- the presence of both dissipations k = k' = 0.05a; 
(Fig.Sa); in the absence of field dissipation k = 0, k' = 0.05a; (Fig.Sb); in the absence of 
molecular dissipation k = 0.05a;, k' = (Fig.Sc). 



Figure 4. The time evolution of Fy{t) as a function of scaled time Lvt for three values 
oiN,k = k' = 0.05a;, a = 1 when 5„,^ / I: = 20, m = 10 (Fig.4a); = 50, m = 25 
(Fig.4b); N = 100, m = 50 (Fig.4c). 



Figure 5. The time evolution of Fy{t) as a function of scaled time a;f for a = 1, when 
Bn,m = 1: N = l (Fig.5a); N = IQ (Fig.5b); N = 100 (Fig.5c). 



Figure 6. The time evolution of Fy{t) as a function of scaled time uot, for a = 1, 
k' = {),k = 0.05a;, when Bn,m 7^ 1: iV = 30, m = 15 (Fig.6a); = 100, m = 50 (Fig.6b). 
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